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Abstract 

Cayley rational forms for rotations are given as explicit matrix polynomials for any quan¬ 
tized spin j. The results are compared to the Curtright-Fairlie-Zachos matrix polynomials for 
rotations represented as exponentials. 

Curtright, Fairlie, and Zachos (CFZ) obtained an elegant and intuitive result |1] expressing 
the exponential form for a rotation about an axis n, for any quantized angular momentum j. as 
a polynomial of order 2 j in the corresponding (2j + 1) x (2 j + 1) spin matrices n • J, thereby 
culminating several decades of previous studies on this or closely related problems m-m- For each 
angle-dependent coefficient of the polynomial, the explicit formula found by CFZ involves nothing 
more complicated than a truncated series expansion for a power of the arcsin function. Although 
a detailed proof of the CFZ result is not exhibited in JT] , the essential ingredients needed to provide 
such a proof are in that paper, and indeed, details of two elementary derivations were subsequently 
given in [7]. 

I discuss here the corresponding polynomial result for the Cayley rational form of an irreducible, 
unitary SU(2) representation with spin j. In this case, the explicit polynomial coefficients involve 
nothing more complicated than truncated series expansions for simple, finite products. I compare 
the Cayley form to the CFZ result. 

The CFZ formula for a rotation through an angle 9 about an axis n, valid for any spin j € 
{0, i,!,§,■■■}, is given by the unitary SU (2) matrix 


2j 1 

U = exp (i 9 h • J) = — A(9) (2 i h • J ) 

k\ 


k =0 


(i) 


where the angle-dependent coefficients of the various spin matrix powers are explicitly given by 


(9) = sin fc (9/2) (cos (9/2))^’^ Trunc 

Lj-fc/2j 


(VI - V 


( arcsin Jx 


( 2 ) 


J z=sin 2 (0/2) 


Here, • J is the integer-valued floor function while e (j. k ) is a binary-valued function of 2 j — k 
that distinguishes even and odd integers: e (j, k ) = 0 for even 2 j — k, and e (j, k) = 1 for odd 
2 j — k. More importantly, Trunc [/ (x)l is the nth-order Taylor polynomial truncation for any 

n 

f ( x ) admitting a power series representation: 


f ( x ) = X f™ xm > 

m =0 


Trunc [/ (x)] = ^ Z™ 1 " 1 • 

m =0 


( 3 ) 
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To emulate the CFZ formula, I now construct Cayley transforms mm of elements in the su ( 2 ) 
algebra as spin matrix polynomials. 

At first glance the Cayley transform for any spin representation would seem to follow imme¬ 
diately from the CFZ result just by changing variables from 9 to a function of the angle, a (9). 
For any given numerical value of h • J this would be so, of course, but it is not obviously so for 
matrix-valued n • J. In fact, it turns out that the matrix produced by the Cayley transform of an 
su (2) generator reduces to a spin matrix polynomial that has a simpler form than the CFZ result 
JT]) for the exponential of that generator. 

For a finite dimensional, unitary, irreducible representation of SU (2), all the underlying spin 
matrices iS = 2 i h • J are anti-hermitian, and for spin j satisfy the 
mm appropriate for (2 j + 1 ) x (2 j + 1 ) matrices. Consequently the 
unitary SU ( 2 ) group element for spin j can be reduced to a spin matrix polynomial [TT] : 


Cayley-Hamilton theorem 
Cayley rational form of a 


1 + 2 ia h • J 
1 — 2 ia h • J 


2j 

E a * ] («) ( 2?: * ■ J ) k 

k =0 


(4) 


where a is a real parameter, and where the coefficients 2 l|^ (a) are to be determined as functions of 
a. The challenge here is to rewrite the geometric series 1/ (1 — 2 ia h • J ) for spin j as a polynomial 
in h • J. Thus define 


1 

1 — 2 ia h • J 


2 j 

E («) ( 2i * • J ) k • 

k =0 


(5) 


Then clearly 2l[j^ = 225^ — 1 and 2ljE = 223^^. The coefficients in the latter expansion follow 
directly from the methods in mm- The results are succinctly given by 


(a) 


det (1 — 2 ia h • J ) 


Trunc [det (1 — 2ia h • J )] , 
2 j-k 


( 6 ) 


where the truncation is in powers of a, and where the determinant for spin j is 


2 j U+1/2J 

det (1 — 2 ia h • J) = (1 — 2ia (j — m )) = fl + 4a 2 (j + 1 — n) 2 J . (7) 

m —0 n=l 

These results are readily checked for small values of j upon using explicit matrices, say h • J = J 3 . 
Indeed, this is how Q was first deduced. 

Subsequently, after the first version of this paper was posted on the arXiv, a detailed proof of 
© was given in m- But, as it turns out, after the previous version of this paper was posted, 
I learned that © and © are true for any N X N matrix M, upon replacing h • J —» M and 
2j + 1 —>• N. The result is a widely known theorem from the resolvent formalism (e.g. see im 
and is proven in numerous textbooks and articles on matrix theory (e.g. see PHI and [15]). 

Still, some comments are in order. Firstly, note that for either bosonic (integer) or fermionic 
(semi-integer) spins, only even powers of a with positive coefficients are produced by the determi¬ 
nant factors in ©. Consequently the 21^ (a) and 23 (a) coefficients have no singularities for real 
a. Secondly, the determinants in ([ 6 ]) are essentially generating functions of the central factorial 
numbers t (m, n ) (see EE]), a fact already exploited in mm- For example, for integer j, 

Li-fc/2j 

Trunc [det (1 — 2ia h • J)] = E 4 m a 2m \t (2j + 2, 2j + 2 — 2m)| , ( 8 ) 
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with the full determinant obtained for k = 0. Thirdly, as j —» oo for any fixed k the truncation in ([ 6 ]) 
is lifted — albeit not without some subtleties m — to obtain, for small a, Hindoo 1 («) ~ a k . 
Nevertheless, in contrast to the periodicized 0-monomials found in [1], the large j behavior here 
does not make the periodicity of rotations manifest. To exhibit periodicity even for finite values of 
j, 6 must be expressed, on a case-by-case basis, as cyclometric functions of a , and then the branch 
structure of those cyclometric functions must be invoked. 

Finally, the Hille-Yosida theorem from resolvent theory provides a one-to-one relation between 
d2D and ([ 6 |). That is, for any hermitian N x N matrix, a Laplace transform gives 


1 

1 — itM 



e~ s e istM ds . 


(9) 


So for linearly independent powers M k , 0 < k < N, as is the case for spin matrices, the matrix 
polynomial expansion coefficients are related order-by-order. 


-j N 1 N 1 ✓»oo 

t _ .TTT = ^2 A n (t) (iM) n , exp ^ C n (</) ( iM) n , A n ( t) = / e~ s C n ( st ) ds . 

1 n= 0 n=0 

(10) 

This provides yet another route to prove the CFZ result (jT|) starting from the known coefficients 
for the resolvent (see m ). For example, for spin j = 1/2, 


exp (i icj) cr 3 ) 



I + i sin (jfj 


03 , 


( 11 ) 


where / 


/ e s cos (st/ 2 ) ds = 

' 0 i w 4 


i r 

1 + I t 2 Jo 


e s sin (st/2) ds = 




1 + 1 t 2 ’ 


det (l - 1 it 03 ) = 1 + 1 t 2 , 


1 


1 


1 _i if(r3 - = TTT75h+ ii + • 

J M and the Pauli matrix is ct 3 = ( J _° x ), as usual. 


( 12 ) 

(13) 

(14) 
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